Electron propagation in crossed magnetic and 
electric fields 



T. Kramerj, C. BracherJ, and M. Kleberj 

t Physik-Departmcnt T30, Technische Universitat Miinchen, 
James-Franck-Strafie, 85747 Garching, Germany 

X Department of Physics and Atmospheric Science, Dalhousie University, 
Halifax, N.S. B3H 3J5, Canada 

E-mail: tkrernierOph.tum.de 

Abstract. Laser— atom interaction can be an efficient mechanism for the 
production of coherent electrons. We analyze the dynamics of monoenergetic 
electrons in the presence of uniform, perpendicular magnetic and electric fields. 
The Green function technique is used to derive analytic results for the field- 
induced quantum mechanical drift motion of i) single electrons and ii) a dilute 
Fermi gas of electrons. The method yields the drift current and, at the same 
time it allows us to quantitatively establish the broadening of the (magnetic) 
Landau levels due to the electric field: Level number k is split into A: + 1 sublevels 
that render the fcth oscillator eigenstate in energy space. Adjacent Landau levels 
will overlap if the electric field exceeds a critical strength. Our observations are 
relevant for quantum Hall configurations whenever electric field effects should be 
taken into account. 



PACS numbers: 73.20.At,03.75.-b,73.43.Cd 
1. Introduction 

The quasi-free propagation of electrons in strong external laser fields is known to play 
an important role for the interpretation of phenomena that cannot be described in 
perturbation theory: The quantum propagation of photoelectrons in a laser field is 
essential for the interpretation of non-linear effects such as the generation of plateaus 
in high-harmonic generation El • More generally, one expects that a study of the 
quantum motion of matter waves in external fields will shed light on the occurrence of 
interesting interference phenomena. In addition, external fields offer a useful testing 
ground for analyzing the interplay between classical and quantum dynamics. 

An impressive example for the application of a static external electric field in 
basic laser-atom physics is the so-called photodetachment microscope 0] which is 
used to determine electron affinities of negative ions with high accuracy. Similarly 
photodetachment studies in a magnetic field reveal oscillations in the photocurrent 
spectrum as a function of the energy of the emitted photoelectrons (for a collection 
of references and new experiments see 0). The observed oscillatory structure is 
clearly related to the existence of Landau levels. In this paper we will analyze the 
quantum motion of electrons in crossed static magnetic and electric fields. One expects 
the subtle properties of a magnetic field to show up most clearly in two-dimensional 
motion in a plane perpendicular to the magnetic field. In this way one generates a 
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Hall configuration with the electric field £ (Hall field) orthogonal to the magnetic field 
B. 

A convenient way of achieving this goal consists in utilizing the method of Green 
functions. The Green function G(r, r';£') is linked to the familiar time-dependent 
quantum propagator K{r,t\r',0) |S] by means of a Laplace transform j7]: 

1 f°° 

G{v,v'-E)^- dte'^*/''if(r,i|r',0), (1) 
i« Jo 

where the propagator represents the probability amplitude for a particle to travel 
from r to r' in a fixed time t. The energy Green function corresponds to the same 
quantum travel, however with fixed energy E. It is important to realize that will 
automatically lead to a retarded Green function, with a matter wave being generated 
at r' with part of it subsequently travelling to r. For free particles the Green function 
is the well-known outgoing wave, 

r'- n - ™ exp(zfc|r-r'|) 



with k — \j2mE /h being the wave number. In external fields, Q will be modified as 
the wave leaves the source point r'. If there are K classical (cZ) trajectories connecting 
r' and r, we have to sum over K amplitudes in order to obtain the Green function in 
semiclassical (sc) approximation: 

K , . s 

G,,(r,r';i?) = ^Afcexp(lM^,'i(r,r';i?)j . (3) 
fc=0 ^ ' 

For the quadratic potentials considered here, the amplitudes A-^ are independent 
of r and r' and are not written down explicitly. What matters here is Hamilton's 
characteristic function W^'5(r, r'; E) for path number fc. W is also known as the reduced 
classical action. It is related to the classical action Sci{v,r' -jtk) via the Legendre 
transformation: 

Sciir, r';tk)^ W^iir, v'-E)- Et^. (4) 

Here tk is the travel time along the classical path number k. Since the Green function 
can be interpreted as the wave function of the moving electron, with r' being the 
starting point, we can readily calculate the current density jsc(r) of the electron. 
According to the laws of classical mechanics, the classical momentum of a particle can 
be derived from Hamilton's characteristic function: 

p^,-VrPFi;(r,r';i?). (5) 

It is well known that cycloids describe the classical motion of classical electrons in a 
combined electric and magnetic field. However, the average classical velocity of an 
electron moving in a plane orthogonal to the S-ficld is the drift velocity, | | — \£\/\B\, 
which is the same for each path k. In view of Q we may therefore write 

{U{v))^e{p{v))xa, (6) 

where j and piv) — \Gsciy, r'; E)^ have been properly averaged. 

What we have demonstrated here is that, by calculating the Green function with 
outgoing-wave boundary conditions we will automatically obtain the £ x S-drift in the 
classical limit. As will be shown below, the Green function also contains information 
to what extent an electron with energy E will participate in the drift motion. To put 
it in a more formal way: By utilizing the appropriate Green function we can calculate 
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the local density of states n{r, E), which plays a central role for calculating the spatial 
flow pattern of matter waves in external fields. To see this more clearly, we rewrite 
the operator equation for G, 

iE-H)Gir,r';E) = S{r-v') (7) 

in terms of a formal solution 

1 



G{v,r'-E)= (r 



E - H 



(8) 



where e = 0^ takes care of the outgoing boundary condition. An eigenstate expansion 
of (m then yields 

where 5';(r) = (r|^/) is the Ith eigenfunction of H with energy Ei. Defining the 
so-called local density of states 

n(r,i?) :=^|vl.,(r)|2<5(i?-i?0 (10) 

we obtain from ^ 

n{v,E) =--lm{Gir,r;E)} (11) 

TT 

It is the imaginary part of the Green function which governs the current of the 
corresponding quantum motion. Therefore, the magnitude of n(r, E) determines to 
which extent an electron that travels through r with an energy E will participate in 
the motion (here the £ x B drift). 

From the gauge property of the Green function in static electromagnetic fields. 



G£,B(r,r';£') = exp 



G£,B{r~r',o-E + qr' ■e){12) 



with g = — e in our case. In the following we will call 

n{E) =n{r = o,E) (13) 

the density of states because n{r,E) is readily obtained from (|12|l . In this paper we 
will calculate the density of states n{E) for two-dimensional motion in perpendicular, 
homogeneous electric and magnetic fields. All necessary results will be derived here; 
they are based on a more general approach that is described in detail in [S] |^ . For a 
wave function approach we would like to refer to 1(1. . 

2. Electric field effects in two dimensional motion 

In the absence of an electric field, the three-dimensional Green function for a uniform 
magnetic field was obtained in closed form in jllj . Also, for a purely magnetic field, 
the two-dimensional density of states is known |12[I13) to show a spike-like structure, 
formally written as a superposition of discrete ^-distributions positioned at the Landau 
levels at E ^ (2k + l)?iWi, where = eB/{2m) denotes the Larmor frequency: 

(^) = ^ E - [2fc + 1]) ■ (14) 

fc=0 
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As we will show below, the presence of an electric field broadens the (5-distribution of 
the kth Landau level into a smooth function of overall Gaussian form that is modulated 
into k + 1 distinct sublevels. 

The desired propagator is governed by the Hamiltonian of a 2D-electron in crossed 
fields, 

hL'b = + l^^^l {^^ + y") -r^-F^- pyXujL + PccVi^L, (15) 

where Fj^ = —e£± denotes the electric force in the x-y-plane. The result is |14| : 

- -2^i^r^{^, [..tcoH.^t) - 1]} . (16) 

To account for the effects of the electron spin, we note that the spin adds a constant 
term zL^gfislS = i-^ghujL to the Hamiltonian and thus merely shifts the energy by 
this amount 22 . The spin dependent densities of states become 

n^AE)^n(^E±^ghu;i}j (17) 

and the total density of states including spin can be mapped back to the one without 
spin: n^i{E) = ni{E) +ni{E). Thus we evaluate n{E) and defer the inclusion of spin 
for the moment. 

For the discussion of the density of states in crossed fields it is useful to transform 
the trigonometric functions in the propagator into a sum over Landau levels. This can 
be done using the identity 



exp[-acoth(z)] ^ ^.-o f ^o) ^^^^ 
sinh(z) ^-^ 



-2z ( k+i 



(18) 

which follows after substituting t — exp(— z) in the generating function of the Laguerre 
polynomials L^"-* (z) |15) . The imaginary part of the Green function therefore reads 

1 °° 
neME) = [C'i^lio, o; E)] = ^ nk,£ME), (19) 

where the partial density of states assigned to the fcth Landau level is given by: 

"^^i f°° J. .-rH''/(4h'')-itE,-/hT (0) /r^t^' 



/OO 



2^2 



^ -^""-[Hk(Ek/r)]^ (20) 



2'=+ifc!7r3/2pr 
Here, the level width parameter 

r = FJ (21) 



is related to the magnetic length I ~ y/h/ (eB), and Ek denotes the effective energy 
shift for the fcth level: 

Ek^E-ry{4tiujL)-{2k + l)nujL- (22) 

In eq. (|20|l . we evaluated the integral over t and expressed the result in terms of 
Hermite polynomials Hk(z) |^. The density of states is identical to a sum over the 
probability densities of the eigenstates of a one-dimensional harmonic oscillator 
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It should be emphasized that the osciUator functions depend here on the energy E 
and not on a space variable. The dimensionless energy variable reads ^ = E^/r — 
Ek/{F±l), with the fcth eigenstate being centered around the fcth Landau level, apart 
from an overall shift r^/(4?ki;i) = y{^±/^)'^ accounting for the kinetic energy of 
the drifting electrons. This result is also consistent with the mapping of the original 
Hamiltonian in crossed fields in eq. (|15|l to the Hamiltonian of a shifted harmonic 
oscillator. Details of the corresponding canonical transformation and its unitary 
representation are presented in Ref. ^^l- The total contribution of the fcth Landau 
level integrated over energy-space is readily available from the normalization of the 
oscillator eigenstates: 

dE nu^sME) = ^ y = (24) 

This result reflects the quantization of each Landau level in a purely magnetic field, 
as given by eq. (|14|) . In the context with fermionic matter waves, as discussed in the 
next section, there is another important quantity: the energy-integrated density of 
states N{Ep), integrated up to the Fermi energy Ep of the system, 

N{Ef) = / dE n{E). (25) 



(A recursive evaluation scheme allows to express N{Ef) as a sum over error functions, 
but we omit the result here.) 

For each Landau level k, the density of states has a Gaussian envelope with 
width r that is split into fc + 1 intervals by the k simple zeroes ^k,i (j = 1, . . . , /c) 
of the polynomial IIk(0. In fig. ^ we plot the resulting density of states for various 
electric field strengths 8±. For small the overlap between adjacent Landau levels 
is negligible, as the density of states drops off exponentially between them. With 
increasing electric field, the Landau levels broaden and finally coalesce. We infer from 
eq. (|23|l that the classical turning point, = ^2k + 1, provides a practical measure 
for the width of the partial density of states nk,£xB{E)- The classically allowed region 
in energy between two Landau levels fc — 1, fc is then given by the ratio: 

half widths of adjacent levels T{^1% + ff) m / /- \ £^ 

^ ^ ■ ^ ^ ~ 1 + ^2fc + 1 — — . 26 

level spacmg 2hu}L \/eh ^ ' » ' 

The overall width of the modulated Landau levels increases with fc^/^. Note that all 
features of the level density of states, including the nodes, scale linearly in width with 
the electric field We conclude that for small magnetic fields, many Landau levels 
overlap, and the density of states becomes smooth. In the limit S ^ 0, an exact 
expression for the two-dimensional density of states in a purely electric field may be 
extracted from eq. (|15|l : 

n£(iJ) = -i3[Gr)(o,o;£;)l 



TT 



27rft2 



Ail (-AE) 



(27) 



where A = 2[m/(Se£_L)2]^/'^, and Aii(z) — dxAi(x) denotes the integral of the 
Airy function |15[ 1^. If additionally £x ^ 0, the constant density of states of a free 
two-dimensional electron gas is recovered: ni,:cc{E) = Q(E) ■ m/{2Tifi^). 
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Figure 1. Two-dimensional density of states n{E) (in units of eB/ (2TTf? uji^)) 
and integrated density of states N{E) (in units of eB/{2-KK)) at four different 
electric fields £^ = 2000, 4000, 8000, f2000 V/m and for a magnetic field H = 5 T 
as a function of the scaled energy E/{fiWL) according to eq. 1201 . Near the fcth 
Landau level at i? = {2k + l)hw]^, the density of states renders the probability 
distribution of a one-dimensional harmonic oscillator in the fcth eigenstate. 



3. Quantum motion of a dilute Fermi gas in two dimensions 

For a non-interacting, i.e. dilute, gas of charged fermions that move against 
a uniformly charged background, the Pauli principle prevents more than single 
occupancy of each energy quantum state. For such a dilute low-temperature gas, 
the available states are filled up to the Fermi level. Hence, the number of occupied 
states is given by (|25f) . Since the electrons arc flowing at a right angle to the electric 
field and not parallel to £±, it is possible to introduce a local Fermi level ^7] 

Epir) = Ef{o) - qr ■ € (28) 

which takes into account that all energy levels are tilted in the direction of the electric 
field. In view of H12|l we conclude that there is translational invariance in the sense 
that 

G£Mr,r;EFir))= GsMo,o:Ef{o)) . (29) 

The exponential suppression of the density of states between Landau levels has 
profound implications on the conductivity in a two-dimensional system, including the 
resistivity plateaus observed in the quantum Hall effect (QHE), as we show now. In a 
simple Drude-like model we may model the dynamics of the electrons by the electric 
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and Lorentz forces amended with a term that incorporates elastic scattering via a 
relaxation time t{E): 

dv TTt 

m— =e£" + evx6 v. (30) 

at T 

Under stationary conditions v is constant and together with the current density 

j — nev the components of the conductivity tensor j = aS become jl3| 

eME)riEl 1 / 1 -.cr^E) \ 

^ ' m 1 + ojIt{EY V ^ct{E) 1 y ' ^ ' 

where loc = — eB/m. For T — > 0, the total conductivity is obtained by 

integrating over the occupied energy range 

dEa{E). (32) 

The availabihty of empty states hmits significant scattering to an energy range of 
order fc^T around Ep. For energies E far from the Fermi energy, t(E) — > oo. In 
strong magnetic fields, the relaxation time t{E) thus satisfies [oJc'TiE)] ^ 1 and the 
transversal component axy mirrors the integrated density of states N{Ep) 

whereas the longitudinal conductivity axx is dominated by contributions in the vicinity 
of the Fermi energy: 



■■'{E)oJct{E) e ujcTjEp) , . 

Kl TTWv2'^'^BT-n{EF)— 2— T^T^- ^4 

BJ-oo ^+^c'^{EY B I + uj^tIEfY 

Thus, the longitudinal conductivity is proportional to the density of states at Ep. 
Finally, we note that the resistivity tensor is related to axx and axy via: 

^ XX ' xy ^ XX ' xy 

From eq. H24(l . it is clear that p^y is quantized between Landau levels. If Ep lies in 
the gap between the {k — l)th and kth level, n{Ep) « is exponentially suppressed, 
N{Ep) = keB/{2TTh) and therefore 

2TTh , . 

Pxx'^O, P^y^l^' (36) 

where, at T 0, the deviations are of the order exp(— /i^w£/r^) = 

exp[— e?i;B'^/(4TO^£^)]. It is interesting to see how single-particle effects can contribute 
to the existence of conductivity quantization. 



3.1. Non-integer filling factors 

We now exploit the fact that the density of states does not only vanish between Landau 
levels, but also within these levels whenever Ep coincides with one of the zeroes ^k.j 
of the individual Hermite polynomials Hk (Ek/F) in eq. (jSOJ- In the vicinity of these 
zeroes, an additional plateau-like structure in the Hall resistivity pxy will appear. 
There, the suppression of the density of states is only quadratic in E. Therefore, 
these plateaus behave differently from the ones due to the integer quantum Hall effect 
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Magnetic field B [Tesla] 

Figure 2. Quantum Hall effect at strong magnetic fields {B > 1 Tesla) for a non- 
interacting two-dimensional electron gas. Hall resistance pxy and longitudinal 
resistance pxx as a function of the magnetic field B for fixed Fermi energy 
{Ep = 0.868 meV). Effective mass m* = 1, effective g-factor g = i, current 
density j^, = 0.2 A/m, t{Ep) = 10~^^ s, T = 0.1 K. The dashed line represents 
the classical Hall resistance pxy with a constant level density. Similar experimental 
results are shown i.e. in 1181 . However, we assumed a high current density in order 
to show the splitting of the Landau levels clearly. 



in eq. (|36|l . Continuing the notation employed above, the new set of plateaus is 
characterized by non-integer filling factors n: 

Pxy = 2 • ^" 

The values of k are obtained by integrating the density of states up to the zeroes of the 
Hermite polynomial. A non-integer k is closely related to the probability of finding 
a particle between the nodes of the fcth oscillator eigenstate in eq. 123|l . Adding 
^fe,o = — oo, S.k,k+i = oo to the list of zeroes of Ufe(^), the contribution Akj of the 
interval between two neighbouring nodes in the kth level becomes 

IwfcCOl^) yielding ^ = fc + ^ A^^i . (38) 

Note that the sum X]i=^i ^k.i — 1 represents the entire level k and thus must attain 
unit value. For j = 0, Kkfi = k, and the integer quantum Hall effect in eq. 
is recovered, while the condition j > yields fractional values for Kkj- The parity 
(—1)'^ of the functions Uk{£,) implies Ufc(O) = for odd k, and thus the existence of 
half-integer filling factors: ki^i — 3/2, K3^2 — 7/2, etc. The other intervals lead to 
irrational filling factors that, however, are surprisingly close to simple fractions. They 
occur for k>2. In the third Landau level (fc = 2), we obtain with ^2,2 — -^(,2,1 — l/v^ 
from eq. ij^ : 



A2,i = A2,3 = ^= + T^erfc ^ = 0.400626 . . . , (39) 




Electron propagation in crossed magnetic and electric fields 



9 



(where e — 2.718281 . . . and erfc(x) denotes an error function JS]). Hence A2,2 = 
1 — 2A2_i = 0.198748 . . .. The sphttings are very well approximated by the rational 
sequence (|, i, |). Similarly we get for /c = 3 A^^i = A^^^ = 0.349992..., and 
A3_2 = As^a = 0.150007 . . . leading to the sequence J,, ^). The plateaus in 
the integrated density of states and transversal resistivity, respectively, due to the 
fractional values of k, are easily identified in all figures. 

3.2. Resistivity plots 

Early experiments to measure the Hall resistivity were performed with fixed magnetic 
field while effectively changing the Fermi energy of the system |19| . A more common 
experimental configuration fixes the Fermi energy and records pxx and p^y as a function 
of B 5H|- The dependence is shown in fig. El In both setups, the transversal current 
vanishes {jy = 0), yielding an implicit expression for the Hall field £y: 



that can be solved for £y under the constraint of a constant longitudinal current 
density j^. 

The physics behind (|40|l is that we allow for charge fluctuations due to the 
presence of external particle reservoirs that feed the electron current. It is interesting 
to see how the changing number N{Ep,B) of mobile carriers that participate in the 
transport generates plateaus in the resistivity. They can be regarded as quantum 
fluctuations about the classical Hall resistance (see fig. 2). Fluctuations in N{Ep) may 
in turn lead to changes in Ep, a possibility which we have not taken into account here. 
This problem is difficult; it cannot be solved within the simple approach presented in 
this paper. 

In practice, interpretation of the plots is impeded by the simultaneous presence 
of two spin populations. Depending on the spin splitting, we expect different fractions 
of the resistivity quantum h/e'^ to occur in the QHE. Let us first consider the case 
of a very small spin splitting (5 — > 0). In this case, the spin dependent (integrated) 
density of states differs from the spin independent quantities merely by a factor of two. 
According to eqs. and dTfjl . the principal plateaus of the Hall resistivity due to a 
completed Landau level are then located at pxy = 2nh/{fe^) with / = 2,4,6,..., 
while the additional plateaus related to the nodes of Uk{x) are characterized by 
/ = 3,^,^,1^,7 etc. The fairly different nature of the two sets of plateaus is 
illustrated in figs. Q] and |21 Once the spin degeneracy is lifted {g ^ 0), due to the 
energy sphtting gtuvL in eq. ifTTIl . the Landau levels observed in the density of states 
n^^l{E) effectively double. Integer QHE plateaus in p^y are then seen for all integer /, 
while the onset of the non-integer effect depends on the value of g. E. g., for g — \ the 
list of expected features comprises plateaus at / = 1, 2, |, 3, |, 4, 5, . . . which 

manifest themselves in fig. El 

3. 3. Breakdown 0/ the quantum Hall effect 

From eq. H26|l one concludes that the disappearance of the plateau between two 
adjacent Landau levels k — 1, k happens at a critical Hall field Sent 



PxyiB.Sy.Ep) j. 



(40) 




(41) 



m V2fc- 1 + V2fc + l' 



Electron propagation in crossed magnetic and electric fields 



10 



Strong experimental evidence for this dependence is reported in Refs. |2UI I21L H ^ . 
With increasing electric field the modulated Landau levels broaden and fill in the 
gaps in the density of states separating the integer QHE. The zeroes of the density 
of states inside the Landau levels are able to withstand stronger electric fields (see 
fig. 0. This is particularly true for plateaus of half-integer index (ki^i = |, ^3^2 = |, 
etc.): Their position only marginally depends on the exact value of the Hall field 
£±, which renders them favourable for experimental detection of the effect. Finally, 
if £± considerably exceeds the critical value £crit, many Landau levels overlap, and 
Pxy ^ 2Trh?B / (emEp) asymptotically becomes a linear function of B, thus recovering 
the classical Hall effect. 

4. Conclusion 

We have analytically calculated how the dynamics of electrons that are injected into 
crossed electric and magnetic fields, depends on the field parameters and on the energy 
of the electrons. The quantum mechanical £ x B drift was analyzed in terms of the 
imaginary part of the appropriate Green function which, in turn, can be related to the 
density of states for those electrons that participate in the drift. It was demonstrated 
how the addition of an electric field will lead to important changes in the density of 
states as compared to those of a purely magnetic field: The discrete delta distributions 
are broadened into smooth functions whose width increases linearly with the applied 
electric field. Their shape coincides with the density distribution of the various 
oscillator eigenstates. Hence, the density of states for the fcth Landau level has k 
zeroes that lead to additional plateau-like structures in N{Ep). For a non-interacting 
Fermi gas of electrons the resulting resistivity plots bear remarkable resemblance to 
actual quantum Hall effect data. Finally, the model predicts a breakdown of the 
integer quantum Hall effect at critical electric fields in accordance with experimental 
observations. 
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